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Abstract 

For many years coherent states have been a useful tool for understanding 
fundamental questions in quantum mechanics. Recently, there has been work 
on developing a consistent way of including constraints into the phase space 
path integral that naturally arises in coherent state quantization. This new 
approach has many advantages over other approaches, including the lack 
of any Gribov problems, the independence of gauge fixing, and the ability 
to handle second-class constraints without any ambiguous determinants. In 
this paper, I use this new approach to study some examples of time repa- 
rameterization invariant systems, which are of special interest in the field of 
quantum gravity. 
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1 Introduction 



The coherent state formulation of the path integral JT] has many advantages over 
a conventional Feynman path integral. Because coherent states describe a minimum 
uncertainty wave packet, there is a natural relation between the classical system and 
underlying quantum system. The coherent state path integral is intrinsically superior 
with regards to a canonical coordinate transform which just amounts to relabeling of 
the states ||. Moreover, for coherent state path integral, it is possible to find a well 
regularized path integral measure (the pinned Weiner measure) 

Recent work 0-0 has included constraints into this formulation. The first ap- 
proach was to consider a semi-classical constraint || . This constraint can be inserted 
by hand into each time step in the construction of the path integral. The result is 
the formal path integral where the action is dependent on the total Hamiltonian. 
Klauder || constructed a projection operator that maps states defined on the full 
phase space onto physical states. The resulting path integral is independent of the 
functional form of the Lagrange multiplier term and hence gauge invariant. Later, 
Klauder and Shabanov generalized this approach to a coordinate-free formulation. 

For constrained system, coherent states offer further advantages. Because the 
path integral is regularized, we are not required to gauge fix to remove the infinite 
volume term that normally appears. The result is an averaging over the gauge orbits. 
Without gauge fixing, there are no potential problems with Gribov ambiguities. Also 
we are not required to eliminate second-class constraints |5j , nor is there the possibility 
of an ambiguous determinant for this case. 

In this paper, we will review coherent state quantization of constrained system and 
compare the results of the author's semi-classical construction of the path integral 
[||] with Klauder's projection operator approach ||. Then we will work through 
details of two time-reparameterization invariant systems. The first example will be 
the single harmonic oscillator. Although this is a very simple example it can give us 
some insight into the details of this formulation. The second example is the double 
harmonic oscillator. 

The double harmonic oscillator is an important example in the study of quantum 
gravity ||. This is a good toy model to help understand the "problem of time."[] 
In essence, one harmonic oscillator can be used as a quantum "clock" to measure 
the other oscillator. In terms of coherent state quantization, the double harmonic 
oscillator show the importance of the geometry on the system Q - [1C]. The geometry 
of the phase space determines the natural kinematical operator in which the system 
should be quantized. For this case, the resulting reduced phase space is spherical so 
the kinematical operators are spin-like operators. Also this system has a potential 
Gribov problem that results in a difference between the ground state energies of the 



1 In a time-reparameterization invariant system, such as quantum gravity, the roll of the local time 
coordinate is difficult to understand. For more about this "problem of time," Rovelli has written a 
series of papers |0. 
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reduced and Dirac quantizations. This Gribov problem results because the gauge 
orbits form a twisted bundle over the constraint surface. 



2 Coherent state quantization 

The ordinary phase space has a natural Heisenberg-Weyl algebra structure that 
comes from the symplectic structure. This operator algebra can be used to construct a 
coherent state representation of the Hilbert space that is labeled by the classical phase 
space coordinates. So, we will begin by considering a set of M pairs of Heisenberg 
operator {Pj, Q k }. These operators obey the standard Heisenberg-Weyl commutation 
relations, 

[P j ,Q k ) = -ih6$ j,k = l,...,M. (2.1) 

The coherent state representation is then a unitary representation of the Heisenberg- 
Weyl group acting on some fiducial vector \rj) chosen from the Hilbert space. 

\p } q ) = e- ifip ' q) e-^e^ p \7]). (2.2) 

In most cases, the fiducial vector is chosen such that the coherent state is "physical 
centered," (r]\Pj\i]) = (i]\Q k \r]) = 0. For this reason, we will choice the fiducial vector 
to be the ground state of the harmonic oscillator, \rj) = |0). This set of states does 
not form an orthonormal basis, as is seen in the overlap function, 

(p, q'\p, q) = exp [\p - p\ 2 + \q - q\ 2 ] + ^ \p' ■ q - p ■ q'}\ . (2.3) 

However, for any choice of the fiducial vector |^), they do admit a resolution of unity, 

In addition, these states form an (over) complete set of states on the Hilbert space. 
We can represent any vector in our Hilbert space as a function of (p, q) by defining the 
function to be ip(p, q) = {p,q\ip). The overlap function , for example eqn.( [2.3| ), is the 
reproducing kernel )C(p', q';p, q) = (p r , q'\p, q) on the Hilbert space. This reproducing 
kernel has the following properties: 

v dpjdq 1 



i>(p',q') = J K(p',q';p,q) i>(p,q)Y[^L, (2.5) 

r N dv' da fj 

K(p",q";p,q) = J JC(p",q";p',q') K(j/,<f;p,q) U (2.6) 
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We can construct a path integral based on this representation of the Hilbert space. 
Unlike the normal path integral which integrates over the configurations space (q), 
the coherent state path integral naturally integrates over the phase space (p, q). To 
construct this path integral, let us start with the Hamiltonian evolution between two 
states \p, q) and \p', q'). The matrix element may be broken in N + 1 time steps. Then 
at each time step we can insert a resolution of unity. 



(p',q'\e-^ T \q,p) = (p' , q'^ 6 ^ le'^^p) 



d(i(pi,qi)\pi,qi)(pi,qi\ 



N 



N 



]J d/i(p n , q n ) Yl (p n+1 , q n+1 \e * H \p n , q n ), 



n=l 



n=0 



\p',q') = \pn+i,Qn+i), 



\p,q) = |Po,?o), 



e=(t'-t)/(N + l). (2.7) 



The measure d/i(p, q) is the same as the measure defined in the resolution of unity 
( |2.4j ). In the limit (e — > 0), if the paths are reguarded as continuous and differentiable, 



then we can formally rewrite the above matrix element (\2.7j ) in form of a path integral 
(see p[ for more details), 



J Dn(jp,q)exp^J (pq - H(p,q))dt 



(2.i 



where the symbol H(p, q) = (p, q\H\p, q). In the stationary phase approximation, this 
then leads to the standard Hamilton's equations of motion. 

Unlike ordinary configuration space path integrals, the phase space path integral 
can be given a natural regularization by inserting an additional term into the path 
integral ||. This is done by changing the measure to a pinned Weiner measure. 
This measure originally arose in the study of Brownian motion. The probability den- 
sity of a particle undergoing Brownian motion is governed by the diffusion equation. 
The fundamental solution of the diffusion equation for a flat metric is a spreading 
Gaussian, 



pit"; t> 



1 



2irv(t" - t') 



exp 



{p"-p'Y+{q" -q 
2v(t" - f) 



A2 



(2.9) 



This solution possesses a semi-group structure with the following product rule: 

pit'"; t') = [ dp" dq" pit 1 "; t")p(t"; t'). (2.10) 
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To construct the Weiner measure, we proceed in a similar fashion as we did with the 
construction of the path integral. We can use the product rule (|2.10| ) repeatedly to 
break the time in N + 1 steps. 



N \ / i x N / N 



(q",p") = {q N+ i,p N +i), (q',p) = (?o,Po), e = (f - t')/(N + 1). (2.11) 
Then, in the continuum limit, we have a formal expression for the Weiner measure, 

dHwip, q) = Me~& /^ +< ? 2 dt Dq Dp (2.12) 

Note that the initial and final points of the paths are fixed (or pinned) on the phase 
space. Writing this in a more general way to include other choices of the metric and 
higher dimensions, the measure is 

2 N 

dfJwfa q) = Ne~^ K^r-) dt J] Dp>DqK (2.13) 

3=1 

The measure in ( |2.ij| ) can now be replaced by the formally well defined pinned Weiner 
measure just by the addition of the extra factor, 

e -^I(wT d \ (2.14) 

In the limit v — > oo, this term formally becomes unity and we are left with our original 
path integral (|2.8). Unless the action is explicitly dependent on the the measure of 
the phase space, this Weiner measure is the only place that the geometry of the phase 
space come into play in the path integral. This geometry (as we will see in our second 
example) determines the natural kinematical operators in which the system should 
be quantized (see 0). 



3 Constrained Coherent State 



In this section, I will review the methods for applying first class constraints 
to coherent state path integrals in general, and then restrict to the case of time- 
reparameterization invariant systems. 

To begin with, let us consider a 2M dimensional phase space labeled by coordi- 
nates (pi,q^) where i,j = 1, . . . M. On this phase space, the constraint surface can 
be defined in terms of a system of N equations 4> a {p,q) = 0, where (2M — N) is 
the dimension of the constraint surface. The evolution on the constraint surface is 
generated by the total Hamiltonian, 
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dF 

H T (p,q) = H(p,q) + X a <f )a (P,q), = {F, H T } . (3.1) 

4>a =0 

For this paper, we are only interested in the dynamics of a system where the 
constraint functions, (f) a = 0, are all first class functions, 

{0a,<M~O, a,b= 1,...,N. (3.2) 

We will assume that the time derivatives will not introduce any new (secondary) 
constraints. Therefore the set {4> a } is complete. This also means that the Hamiltonian 
is also a first class function, 

{0 a , H T } = {<f) a , H} + A b {0 a , <M « 0. (3.3) 



dt 

Furthermore, because these commutators are all weakly vanishing, near the constraint 



surface, they are given as linear combinations of the constraint functions [12]. Thus 
the Hamiltonian and the constraints form a closed algebra, 



{H,4> a } = h c a cf> c , (3.4) 
{0 a ,<M = C Q y> c . (3.5) 

We can now go on to study solutions to the time evolution equation on our con- 
straint surface, 

^ = {f,H T } = {f,H} + \*{f,<l> a }. (3.6) 



We see that in general the solution to eqn. fl3.6|) will depend on the choice of the 
Lagrange multiplier A a . However, a physical observable will not have any dependency 
on this choice. Therefore, any solution that differs only by changing the value of the 
Lagrange multiplier is defined to be equivalent. For first class constraints, these gauge 



transformation are generated by the constraint equations [12 



5f = 6e a {f, <f> a }, (3.7) 

and the dimension of the gauge transformations is the same as the number of con- 
straints. Thus the reduced phase space (the manifold after applying the constraints 
and quotienting out the gauge orbits) is then a (2M — 2N) dimensional manifold. 
This reduced phase space admits a local symplectic structure |fT2| . So, it is possible 
to locally find a canonical coordinate system such that 

= 1 V. (3.8) 
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Let us consider, these two phase space. In general we may describe the dynamics 
of the system on either the full phase space in terms of the total Hamiltonian or on 
the reduced phase space in terms of the reduced Hamiltonian, 

H = H(p,q)\ (3.9) 

<Pa=0 

So we already have the a set of coherent states on the the full phase space, 

|p, q ) = e - ifM e-^e* q - p \ri), (3.10) 

where P, Q obeys the standard Heisenberg-Weyl commutation relations [Pj,Q k ] = 
ihSj. So now, let us take the naive approach that we can construct a coherent state 
on the reduced phase space in the following way. We can try to use the symplec- 
tic structure on the reduced phase space to define an other set of Heisenberg-Weyl 
operators [P-', Q k '] = iHS^ where j,k = 1, . . . , N. Note that these operators may 
not be global well defined nor are they necessarily defined in terms of the Heisenberg 
operators from the full phase space. For a least for some covering space of a large 
patch, we can construct the coherent state, 

\p,q) = e-ifode-fr&ei^lri). (3.11) 

If the initial and final states, \p,' , <f) and \p", q"), can be lifted back up onto the full 
phase space, \p',q') and \p",q"), the the resulting dynamics should be equivalent (up 
to possible normalizations), 

(p\q'\e-^'-V\p,q) ~ (^| e -**°<^>|p,g>. (3.12) 

As we can see, there are two problems that we need to deal with in comparing 
these two descriptions. One is that we need to construct a set of meaning coherent 
states on the reduced phase space. We will see how this is done when we consider 
the projection approach of Klauder's ||. The other problem is to understand the 
dynamics of the total Hamiltonian in terms of a system of coherent states. 

To begin with, let us consider the evolution generated by the total Hamiltonian 
(|3.3|) . We can use the resolution of unity on the full phase space to construct the 
path integral ( |2.7| ). This gives us the path integral, 

M J DpDqexp-^y ih(p, q\^\p, q) + (p, q\H + \ a $ a \p, g)|. (3.13) 

We can replace the operators in terms of either the "upper" or "lower" symbol depend- 
ing on our construction of the path integral (see U for more about these symbols). 
The resulting path integral is then 

M J DpDqexp-^ [Jpq- H(p,q) - A a a (p,g)} . (3.14) 
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We see that the results of the path integral depends on the choice of the Lagrange 
multipliers \ a (t). To fix this, we can extend our phase space to include the Lagrange 
multipliers (see for an example of such an extendition) . or reconstruct the path 
integral by placing the constraints in at each time step by hand 0. The resulting 
new path integral includes integrating over the Lagrange multiplier, 

M J DpDqDXexp |~ J pq - H(p,q) - A a a (p,g)| , (3.15) 
and when we integrate over A a , the result is our constraint equation, 

J D\exp{\ a <j )a } = 5(<p a ). (3.16) 

Another method of imposing the constraints is to project the coherent states onto 
the physical states ||. This projection operator for first class constraints commutes 
with the time evolution. Thus a physical state will evolve into another physical state. 
When this projection operator is included, the resulting path integral picks up an 
additional term which is just a normal integration over the Lagrange multiplier. The 
resulting path integral becomes independent of the functional form of the Lagrange 
multiplier. 

In either case, we are still left to contend with the gauge degrees of freedom. 
Normally, when we integrate over these degrees of freedom, we will get the volume 
of the space of paths for the gauge orbits. In an ordinary path integral, this volume 
term would be infinite and we would have to include a gauge fixing term to remove 
this infinite redundancy. With a coherent state path integral, we can use the Weiner 
measure to regularize the path integral, and because of this, we are not forced to 
introduce any gauge fixing into the system. The result is just an well defined averaging 
over the gauge degree of freedom. Then because we not required to gauge fix the 
system, we avoid any possible Gribov problem.^] 

Let us begin a more detailed construction of this path integral by considering 
Klauder's projection operator approach ||. We wish to find a projection operator 
that takes any state onto a state that is annihilated by the constraint operator (or 
physical states), 

\p,q) P hys = P\p,q)- (3-17) 
As a standard projection operator P must have the following properties: 

P 2 = P and F f = P. (3.18) 

We can construct an example of such a projection operator in terms of the constraint 
functions. As we have seen, the constraint functions form a Lie algebra ( |3.4| ). Let 

2 A Gribov problem or obstruction occurs when the the gauge fixing term can not be defined 
globally or that gauge fixing function intersects a gauge orbit more then once. 
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us assume that the we can find the corresponding constraint operators such that this 
algebra is carried over to the commutator algebra. We can use the group elements 
generated by these operators to form a projection operator, 



P= I e iXa ^dfi(\). (3.19) 



We will choose the measure to be normalized / dfi(X) = 1. In addition, it must satisfy 
the above properties of the projection operator fl3.18|) . For a compact group, such a 
measure is the left and right equivalent Haar measure (see || for more details). This 
projection operator then projects onto the states that obey the quantum operator 
equation & a \i>) = 0, which are the physical states. 

For a non-compact group, finding a measure that is normalizable is a bit more 
difficult. Klauder || suggested the following idea. Let the measure take the form, 

F = f e »*(™¥*d\). (3.20) 



VTA 

This projects onto states where the constraints operator is within a small interval, 

\\Mp,q) P hvs\\ < \\e \p,q) P hys\\- (3.21) 

Then in the limit e — > 0, we have a handle on how to regularize this measure. 

In order to use this projection operator in our construction of a path integral, we 
note that because the measure is left invariant, the projection operator is invariant 
under gauge transformations that are generated by the constraint, 

= J e ixAa dn(o- 1 ■ A) 

= P. (3.22) 

In a similar fashion because the addition of the Hamiltonian operator into the algebra 
is also closed (see eqn. |3~^) , the projection operator commutes with the time evolution 
operator, 

Pe~% A = e~% 6 P = Pe~% 6 P. (3.23) 

For the time evolution of the physical states, because the projection operator com- 
mutes with the evolution operator and the Lagrange multiplier term can be absorbed 
into the projection operator, the matrix element can written in terms of just the 
evolution of the physical state on the full phase space. 
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(p, q; t\p, q) phys = (p , q'\e * Ht P\ P , q) 

= (p',q'\e-% 6+iaa * a P\p,q) 

= (p',q'\e-^P\p,q) (3.24) 

Then we can place a resolution of unity between the projection and time evolution 
operator. After doing this, the first term becomes the evolution on the full phase 
space, for which we have already constructed the path integral ( |2.7| ). Formally, we 
have the following modified path integral (see f| for more details): 

J exp^Jpq-H{p,q)^{p",q"\P\p,q)DpDq. (3.25) 

By not absorbing the Lagrange multiplier term into the projection operator, we can 
repeat the same process with the Hamiltonian replaced by the total Hamiltonian 
operator. The resulting path integral is 

Jexp^Jpq- H(p, q) - \y a (p, q)\{p", q"\P\p, q)DpDq. (3.26) 

So even though this path integral appears to still depend on the choice of the func- 
tional form Lagrange multiplier, we see that in fact it is equivalent to the path integral 
without this term (|3.25|) . 

Now let us consider a different construction of the path integral. Let us work with 
only one constraint <f)(p, q) = 0. We have the projection operator given by 

P = J e~ iT$ dp'(T). (3.27) 

Let r(ti,t 2 ) = ft? <^(t)dt. The operators are inherently time independent. So, we 
have 



P 



Je~ l ^ m * d %'(T( tl ,t 2 )). (3.28) 



Then using the properties of the projection operator ( p.!8|) , we can construct the 
simple product rule, 



P= P 2 

= J fe^ A(t) *^^ AW *V(r(ti,fe))dA* / (r(* 2) *3)) 
= | e -r i 3 ^V(r(t 1 ,t 2 ) + r(t 2 ,t 3 )) 
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Jt 1 



rf/i'(r(ti,t 3 ; 



(3.29) 



This projection operator can be broken into N time segments as we did in the con- 
struction of the path integral (p.7|) . So we can repeat this construction to include 
the projection operator. The above constrained propagator ( 3.24 ) then can be writ- 
ten in terms of the discrete path integral. Let each time step be given by e, then 
jtn+ \^(if ~ e \ n anc [ the measure is <i/i'(r) = dp,'(eX). The discrete path integral is 




J] d^{p n ,q n ) J] dti'{e\ n ) ]J (p n+1 , q n+ i\e-% V A "*|p n , q n ). (3.30) 



Let us rescale the Lagrange multiplier A — > X/h. Then we see that have derived the 
time evolution operator in terms of the total Hamiltonian, 



N N / _ \ N 



f f[ MPn,qn) ft d V' (l A n) II>n + l,gn + l|e-^ +A "*V,g n >. (3.31) 
J n=l n=l ^ ' n=0 



In the Continuum limit, we wish to replace the ordinary measure above with the well 
defined Weiner measure. Certainly, we already know to do this for the first measure of 
the momentum and position, but we would also like to do the same for the Lagrange 
multiplier measure. 



Looking carefully at (|3.29|) , we see that in fact we already have a path integral. 
We can also see that if we let the measure be defined in terms of the fundamental 
solution of the diffusion equation (|2.9|), then we can write the measure as 




p(t 1 ,t 2 ) dX(t 2 



1 MAfal-AWq (3.32) 

We see that this measure is normalized, 

J dfi(r) = J p{t u h)dX(t 2 ) = 1, (3.33) 

and the product rule of the this measure ( |2.10| ) is consistent with the above product 
rule Q3.29 ). Then in the formal limit, we should replace the measure with a Weiner 



measure. Note however that this Weiner measure is not pinned at both ends, but in 
fact we should integrate over the end terms. This integration is how the propagator 
loses its dependence on the Lagrange multiplier. So this measure can be taken into 
an unpinned Weiner measure. Formally we can write the path integral as 
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AfJd^(p,q)d^(X)expUj (pq-H(p,q)-\(t)<j>(p,q)dt)}. (3.34) 

This path integral is discussed in earlier work by the author |4| in terms of the semi- 
classical construction of the path integral. 

Now that we have constructed the various forms of the path integral for the con- 
strained system (eqs. |3.25| , |3.34j ), we would like to consider time-reparameterization 



invariant systems. A large class of time-reparameterization invariant systems may 
be written in terms of a single constraint, Ht = X(H — E). A rescaling of the time 
coordinate can be aborted into the definition of the Lagrange multiplier. So X(t) is 
just a lapse function. 

In terms of the coherent state quantization, the matrix element is quite simple. 
Because the Hamiltonian is zero there is no "time" evolution on the full phase space, 
the matrix element (|3.24j) is then just 



(p', q'\ t\p, q) phys = (p'q'\P\p, q) 



AH-E) 



d^r){p",q"\e-^^\p',q'). (3.35) 

Because there is no dependence on the position and momentum through the Hamil- 
tonian, the path integral ( |3.34p becomes trival to integrate in this direction. The 
remaining path integral is just dependent on the Lagrange multiplier. Then we can 
use the product rule ( |3.29| ) to integrate along this direction. The resulting matrix 



element is the same as above ( |3.35|) . Note that the integration variable for this oper- 



ator should be identified with the proper time t = Hi Ht)dt This was first noted by 
Govaerts |J in his consideration of the free particle case. 

Let us now consider two examples of time-reparameterization invariant systems; 
the single and double harmonic oscillator. 



4 The single harmonic oscillator 

The harmonic oscillator is a natural place to begin the study of time reparame- 
terization invariant systems. In addition to being a simple system to work with it is 
also the natural setting in which coherent states first appeared.^ In this section, we 
will compare the recent projection operator approach to standard Dirac and reduced 
phase space quantization. 

We will begin with a quick review of the classical time reparameterization invariant 
harmonic oscillator. The total Hamiltonian for this model is given by 

3 In 1926, Schrodinger was interested in finding a wavefunction for the harmonic oscillator where 
the center oscillated at the classical frequency [fL4|. 
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Ht — X 

The action for this system is then 



l(p 2 + u;y)-4 (4.1) 



pdq- J X (p 2 + uo 2 q 2 ) -E~\dt. (4.2) 

The equations of motion can easily be solved by defining the proper time r = Jq Xdt. 
Then the equations of motion appear as the normal equations of motion for the 
harmonic oscillator with r replacing the time variable 

d P 2 d q u o \ 

Tr=- Uq Tr =P - (4 - 3) 
The solutions of these equations of motion are 

q = A cos(co>r + cf)) p = Au sm(u)T + <p) . (4-4) 
In addition the equations of motion we also must satisfy the constraint equation, 

^(p 2 + u 2 q 2 ) -E = 0. (4.5) 



Substituting the equations of motion ( [4.4[ ) into the constraint equation, we can solve 
for the amplitude, 



A = (4.6) 

The remaining degree of freedom of this system is just the gauge degree of 
freedom. To see this, let A — » A + e, then (f> — » <// = (ft + et + 0(e 2 ). So the resulting 
reduced phase space is just a single point. Quantizing this system is trivial, because 
there is only one state. Note, however, that the energy E appears to be arbitrary. 

Now, let us look at the Dirac quantization of this system. To begin with, we must 
find the operator corresponding to the constraint function. It is natural to choose 
a Hermitian operator. For convenience, we will switch to the complex coordinate 
a = sf^hQ + *\/ 2uhP> anc ^ we w ^ re pl ace the momentum and position operators by 
the standard harmonic oscillator raising and lower operator (a,a^). The constraint 
operator can be written 

$ = ^( aa t + a )a)-El. (4.7) 

Let us define E' = E/uTi — 1/2 and rescale A. Then the constraint operator can be 
written in terms of the number operator (era), 

$ = a f a - E'l. (4.8) 
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Following Dirac quantization, the physical states are defined as the states that are 
annihilated by the constraint operator. Therefore, the physical state is an eigenstate 
of the number operator. This also imposes the condition that E' is an integer. 

a f a|*W = #'I*W l*W = l™>> E ' = n - ( 4 - 9 ) 

The Dirac quantization also leads to the single state \n). However it imposes the 
restriction that the energy is quantized E = hu(n + 1/2). 

We would now like to consider this system in terms of coherent states. We can 
project the coherent state on the full phase space \a) onto the physical space by using 
the projection operator P. 



Ot)phys = 



„ / oo n 

' V n=0 Vn\ 



n) 



oo 

l«| 2 /2 



£ ^ (/ e^ n - E '^\^ \n) (4.10) 
Using the measure for non-compact groups ( |3.20|) , the integral becomes 



e 

n=0 



/ 







E' — 


n 


dX = < 


} 


E' — 


n 






E' - 


n 



e iX(n-E') (^^-) dX = {1/2 \E'-n\=e. (4.11) 



We can choose e to be arbitrarily small. Therefore, we see that E' must be arbitrarily 
close to an integer m otherwise the physical vector is null. So, if we let E' = m, we 
can calculate the above sum ( |4.10| ). 



|aW = e- |Q|2/2 ^=|m> (4.12) 
ym! 

However, this physical state is not yet normalized in the new space. After normalizing, 
the physical state is the energy eigenstate with a phase factor out in front. 

\ a )' = J^Vr = rrH™> = eMn V>- (4-13) 

1 /phys \(a\a)\ phys \a\ ml 

This phase factor is obviously irrelevant to the physics of this system, and it is easy 
to see that, in fact, it is the gauge degree of freedom generated by the constraint 



It is clear that the projection method in this system is equivalent to the Dirac 



quantization ( fl9l) . 



The "evolution" of this system in the reduced phase space is trivial since there is 
only one state. However we would like consider the matrix element on the full phase 
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space so we can compare with the classical solutions. On the full phase space, we can 
look at the physical state on the full space space. Let the physical state be labeled 
by 



\<p a >) = \a f ) phys , (4.14) 

where \a') is the state that is normalized in terms of the full phase space (|4.12|) . The 
matrix element, which is also the wavefunction on the full phase space, is then 



<M«") = (a"\<f> a ,) = (a"\e-\ a '\ 2 / 2 ^=\m) 



/m 

a 


\m) 






'| 2 /2( 


a'a") m 


m! 



(4.15) 



Let a' = r'e and a" = r"e , and let us renormalize the wavefunction such that it 



is approximately one at the peak, r 



/2 r/2 
—- r = m. 



y'[Oi 



V2 



tm t/m 
-|r"| 2 /2 -|r'| 2 /2 r r 

Tim e ' e 1 1 ' : — . 

ml 



(4.16) 



This normalization can be explained in term of a "gauge fix". We want the phase 
factor from the initial physical state (|4.13 ) to be one e md = 1. Then, normalizing 
this function ( |4.16j ) over phase space, we have 



<f>a'(a")\ 2 \5 ( 



Amd' 



1 



' da'da'\ ( da" da 1 ' 



71 



7T 



(4.17) 



So there seems to be a natural choice for the form of the gauge fixing term in this 
system. 

Ordinarily, we would have to construct a set of gauge invariant operators to work 
on the reduced phase space, but in this system, the gauge orbits are are understood 
to be the phase of the state. Because we know the behavior of the gauge orbits, we 
can remain on the full phase space and consider the correlation between the physical 
state and the other coherent states on this space. Some of the important correlation 
functions of this system are 



(a"\H\<t> a/ ) = (™ + £) (a?\<f>a>) 

(«"|Q|^) = \/^(^ + «") <«"K') 
(a"|P|0 Q ,) = iV2?m (^j - a"^j (a"\(f> al ). 



(4.18) 

(4.19) 
(4.20) 



14 



Once again these correlations are peeked when the classical classical constraint func- 
tions are meet, \a"\ = \a'\ = m. Then, the classical limit (expanding about the peak) 
gives 



(a"\H\<t> a .) = m + 1/2 + 0(h), (4.21) 

(a"\Q\<f>a>) = — cos(6" - 9') + 0(h), (4.22) 

(a"\P\(j) a t) = VE sin(0" - 9') + 0(h). (4.23) 

Then, we can identify 9" — 9' = uot + <p. The resulting corrections then give us back 
the classical equations of motion (eqs. [Oj, [PI ). 



5 Double harmonic oscillators 



Next, we would like to consider a system of two independent but identical har- 
monic oscillators. In addition to being a non-trivial example of a time- reparameter- 
ization invariant system (it still has two degrees of freedom remaining after applying 
the constraints), the double harmonic oscillator has been of interest in helping to 
understand the "problem of time" in quantum gravity. One of the oscillators can be 
thought of as a quantum clock. Then the other oscillator can be written in terms of 
the "time" that this clock reads (see || for more about this system). 

In this system, we also encounter a potential Gribov problem. The constraint sur- 
face is topologically a three sphere S 3 . The gauge orbits are topologically equivalent 
to a circle S 1 . The resulting reduced phase space is the two sphere S 2 . However the 
three sphere is not a trivial bundle over the two sphere, S 3 ^ S 2 x S 1 , but rather a 



twisted bundle. Therefore, we can't find a global gauge fixing condition |12|| . We can 
find a local gauge fixing and extend it to cover all but a single point of the gauge 
orbit. How we treat this point will determine the ground state energy for the reduced 
phase space quantization. 

Once again, let us start by considering the classical system. We will choose each 
of the harmonic oscillators to have the same frequency uj\ = 002 = to. Then the 
Hamiltonian for the double harmonic oscillator is given by 

Ht = X [\(Pi 2 + wV) + \( V 2 2 + ^V) - • (5.1) 
The action of this system is given by 

S = Pxdqi +p 2 dq 2 — H T dt. (5.2) 



As in the single harmonic oscillator case ( f4.3|) , the equations of motion are easily 
solved in terms of the proper time r = f X(t)dt . 
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qi = A cos(ur + 0) Pi = Auj sin(u;r + 0) 

q 2 = B cos(c<jr + 0') p 2 = Bui sin(ur + 0'). (5.3) 



The constraint equation, 



^( Pl 2 + o/V) + ^(P2 2 + w V) = E, (5.4) 



limits the amplitudes to 



{Auif + (Bui) 2 = E. (5.5) 



If A — > A+e, the we see that the gauge transformation take both — > (j)+uiet and 0' — > 
0' + wet. So the degree of freedom that is independent of the gauge transformation is 
the difference between the initial phase of the two harmonic oscillator A0 = — 0'. 
The resulting reduced phase space is two dimensional. 

We can take for the coordinates on the reduced phase the momentum and posi- 
tion of the first harmonic oscillator. This set of coordinates inherits the symplectic 
structure from the full phase space {qi,pi} = 1. However, the metric on this reduced 
phase space is no longer flat. Let us rescale the the momentum pi — > \fuSh p, t and the 



position qi — * J%/ui qi such that they have the same units. Then the volume of the 
full phase space is given by 

™=/n(^V (5..) 



i=i 



We have the standard Cartesian metric, 



da 2 = dqi + dpi + dq 2 + dp 2 . (5.7) 

To find the metric of the reduced phase space, we will restrict the coordinates 
to the constraint equation and choose the local gauge fixing term (tan -1 (p 2 /<?2) = 
constant). This gauge fixing is not global because it is ill defined at p 2 = q 2 = 0. In 
order to find the induced metric on the reduced phase space, it is easier to work in a 
set of two polar coordinates, 



2 2,2 2 2,2 

n =p\ +qi r 2 =p 2 +q 2 

6 X = tan" 1 ( ^ J 9 2 = tan' 1 (f) (5i 



The constraint equation Q5.4D , after our rescaling and change of coordinates, looks 
like 
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S 8 = n 2 + r 2 2 = (5.9) 
urn 

After apply the constraint (rxdri = —r 2 dr 2 ) and the gauge fixing (d# 2 = 0), the 
metric on the reduced phase space becomes 

da' 2 = (l - ^pj dn 2 + r 1 2 d6 1 2 . (5.10) 

We see that this metric is a constant curvature metric (R = 2/S 2 ). The metric is 
ill defined at r\ = S, which is also the same place that the gauge fixing term is ill 
defined (r 2 = p 2 = q 2 = 0). 

We can use this induced metric on the reduced phase space to tell us a bit more 
about the system. Let us follow a similar system that Klauder discussed [||-Q On the 
two sphere the total surface area must be quantized in order that the term exp(i § pdq) 
for a closed path be unambiguous. Note, we can not include the possibility that S = 
because the metric (and the gauge fixing) are ill defined. Hence 



27m 



J dp A dq = J y/g dpdq = irS 2 n = 1, 2, 3, . . . (5.1 1) 

This implies that the energy is also quantized E = hum where n = 1,2, 3,.... The 
reduced Hamiltonian is just zero so the resulting propagator on the reduced phase 
space depends only on the Weiner measure. We can see this in the path integral 

M J Dp 1 Dq 1 e-&H& 2 e i I pi * tl . (5.12) 

Such a Weiner measure gives rise to spin-like kinematical operators S{ where [Si, Sj] = 
itijkSk |J • We will not work through the details of the resulting spin system here 
because the details can be seen when we consider the Dirac quantization. 

Now we would like to consider the Dirac quantization of this system (see [BJ for 



a similar discussion). Let us begin by constructing the projection operator (|3.18| ). 
Similar to the single harmonic oscillator ( |4.8| ) the constraint operator can be defined 
in terms of the raising and lowering operators for the independent oscillators. 

§ = a ] a + b ] b- E\ E' = E/uh — 1. (5.13) 

Because each oscillator is independent (before the constraint is applied), the double 
harmonics oscillator has a complete set of vectors that is just the direct product of 
the eigenvalue of each of these number operators, 



4 Note, the form of this metric is slightly deferent then Klauder's. The results is a difference of a 
factor of two in terms of the area. 
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\m, n) = \m) <g> \n). (5-14) 

a' a\m, n) = m\m, n) , 6'6|m, n) = n\m, n). (5.15) 

The constraint on this basis then quantizes the energy E = hu(m + n + l). In terms of 
the Dirac quantization the energy would have to be an integer (E = m! = 1, 2, 3, . . .). 
The resulting physical states would be given by 



\^) P hys = \n,m' - n). (5.16) 

On this set of states, the raising and lowering operators form above (5T4j ) are not 
defined, in that (a\<&) p hy S ) is not a physical state. So, we need to find another set of 
operators that are defined on this set of states. These operators are equivalent to the 
spin operators. 



For convenience, let us continue to use the scaled momentum and position (5.6). 
In term of these coordinates, we can define a new set of coordinates which have zero 
Poisson brackets with the constraint {s«, 0} = 0. They are 

«i = ^(PiP2 + qiq2), (5.17) 
S2 = ^{P2qi -Pife), (5.18) 

s 3 = \(pi 2 + qi 2 -P2 2 -q2 2 ). (5.19) 

This set of coordinates possess the standard 50(3) Lie algebra, {si, Sj} = eijkSk- The 
square of these three coordinates is the constraint surface radius 

S! 2 + s 2 2 + s 3 2 = s 2 = ^S 2 . (5.20) 

Because the coordinates have a zero Poisson bracket with the constraint, this set 
of coordinates is gauge invariant. However, they are not all linearly independent, 
so let the reduced phase space be described by si,s 2 . The induced metric is flat 



da' = so 2 (dsi + ds2 ), and the domain is just a disk (si 2 + s 2 2 < s 

We can write a set of operators that correspond to the classical coordinates above, 
eqs. ( |5.17| ) - ( |5.20| ), that preserves the £0(3) algebra in terms terms of our raising 
and lowering operators fl5.14|) . 



Si = -(a6 f + a f &) S 2 = l -{atf -a)b) 

S 3 =-{a t a-tfb) 5 = -(o f o + 6 f 6) (5.21) 

_ 
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Then, it is possible to map the physical states give above ( |5.16| ) onto the set of angular 
momentum eigenstates. 



Sq\u, m! — n) — —m'(m' + l)\n, m' — n) = j(j + l)|j, m), 



Ss\n, m — n) = — m')\n, m! — n) = m\j, m). 



(5.22) 
(5.23) 



This means that j = 2m' and m = n — j . The raising and lowering operators for the 
angular momentum S± = Si ± 1S2 act on this set of states in the normal way, 



S+\j,m) = s] (j -m)(j + m + l)|j, m + 1) 
S-\j,m) = J (j + m)(j-m + m - 1). 



(5.24) 
(5.25) 

We can now construct the SO (3) coherent states form these operators ( p.24| ) (see 
[0 for the details about this coherent state). The coherent state is then 



|£> = exp(6S' + -6S'_)|7 7 >. (5.26) 
The let us choose the lowest weight vector from above (|5.24 ) as our fiducial vector 



1 77) = \j, — j). Then we can rewrite the above coherent state representation (|5.26|) in 
terms of the above basis vectors ( |5.22|) . 



io = (i + i£i 2 r E 



2j! 



,„■—, » (./ • i'i)Uj - 1,1)1 
This coherent state has a resolution of unity, 



e +m \j,m)- 



(5.27) 



(2j + 1) dCd[ 

7T (1+1^ 



lexei- 



(5.2J 



With this resolution of unity, we can construct the path integral. In the continuum 
limit, this path integral appears as 



Af / djj! v w exp < ih 



J 



(5.29) 



The Weiner measure for this system was described in We see that in fact the 
reduced phase and the Dirac quantization result in a spin system, where the energy 
of the system is mapped onto the total angular momentum. 

On first appearances would seem that the reduced phase space and Dirac quantiza- 
tion lead to the same results. In both cases the the energy is quantized E = uh(n + l) 



where n = 0, 1, 2, . . . (see eqs. |5.11| , |5.16| ). However, this result is dependent on the 
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fact that we removed the point where the gauge fixing went bad from the reduced 
phase space. It is possible to include this point by using two coordinate patches in- 
stead of one, where each coordinate patch has its own gauge fixing. Then, we can map 
this gauge fixing across the boundary. In so doing, we can include the zero energy 
term in the reduced phase space (|5.11|) . The result is that the two systems then have 
different ground state energies |15| . 

In the case of constrained coherent state path integral, we will end up integrating 
over the gauge orbits in effect averaging over all possible gauge orbits. Because of 
this, we will not have to fix a gauge and we will not encounter this Gribov problems. 

Let us work through the projection operator approach to the constrained coherent 
states for this system. Extending the single state oscillator, the coherent state for the 
double harmonic oscillator can be written as 



\a, (3) = e -H 2 /2-l/3|72 £ / _ a m p n \m, n). (5.30) 

m,n 

Then we can project this on the physical states. 



r / 00 1 

l«,/3W= / e^d^(X) e -N 2 /2-|^/2 £ a m ^\m,n) 
J \ — „ x/nWmA 



m ,n vn!vm! 

00 1 

3 -M 2 /2-|/3| 2 /2' 



< " - ' f] 1 ^ a m (3 n ( [ e iX{n+m - E) dfi(X)) |m, n) (5.31) 
m r, vn! vm\ \J J 



We will again choose Klauder's measure for non-compact groups ( |3.20|) for the mea- 



sure for this projection. Then similar to the single harmonic oscillator (|4.11 ), the 



physical vector is null unless E is arbitrarily close to an integer. So let E = m' 
m + n. Then the physical vector is given by 



/ 1 

\a, f3) phys = e-H 2 -^ 2 ^ ^ — — a ^-"|n, m' - n) (5.32) 

Now we wish to normalize the physical vector. It is just a quick calculation to show 
that the normalized physical vector is 



\a,P) phys = (\a\ 2 + \(3\ 2 ) 2 Yl \ hr T 1 ,, u n P m '- n \n,m'-n) (5.33) 



n=0 



It is easy to see that the gauge transformation generated by the constraint is 
a — > ae t8 and (3 — ► f3e ld . Like the single harmonic oscillator (|4.13|) this gauge trans- 
formation appears as an overall phase in front of the physical vector 

\a, P) ph ys -> e im ' \a, (3) phys (5.34) 
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To remove the gauge dependence, let us then define £ = a//3. Writing the physical 
state in terms of this variable, we can factor out the gauge transformations which 
appear again as phase factor ( |4.13| ). 



\a,0) 



phys 



jm'e 



a 
P 

m' m' 

2 E 

n=0 



E 



m'l 



a 



- i' i/ / i i "n |n,m'-n> 



m'! 



nlfm' — n) 



£ n \n, m' — n). 



(5.35) 



It is easy to see that the physical coherent state maps onto the 5*0(3) coherent state 
(|5.27 ). The energy is mapped onto total angular momentum j = IE' . 



10 = (i + KIT 3 ' E 



2j< 



m=-j 



(j + m)!(j — m)! 



e +m I J, m> 



(5.36) 



The resulting reduced phase space agrees with the information that we were able 
to discern from the reduced phase space coherent state discussed earlier ( |5.12|) . It 
is clearly a spin system with the total angular momentum given by the the energy 
(2j = E>). 

The propagator for this system is simply the overlap function of the 5*0(3) coher- 
ent state, 



(t%) = (a',(3'\P\a[3) 

= (i + iri 2 )- j (i + i£i 2 )^(i+ro 2j - (5.37) 

At any given "time", the coherent state gives a minimum uncertainty wave packet 
(Ji 2 )(J2 2 ) = ^(^o) 2 - The most probable matrix element (the classical solution) is 
simply £' = £. Any expectation value of the spin operators ( |5.21| ) is simple enough 
to calculate as well. From these expectation value the classical "dynamics" of the 
reduced phase space can be can deduced. 

Returning to the full phase space, we would like to reconstruct the the classical 
equations of motion as we did in the single oscillator case ( EL"2l) . Using the non- 



normalized physical vector (|5.32j ), we can define the physical state on the full phase 
space as 

=P\a',(3') (5.38) 
The wave function of this state on the phase space is again the overlap function, 
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<j>a!pi(0i",f}") = (a",/3"|0 a / )( g/) 

= e -|(l«'T+l/3"l 2 +KI 2 +l/3'l 2 ) ft"®' ± 0"0') m / 5 _ 39 n 

ml 

We can renormalize the state such that at its peek, it is approximately one, as we did 
in the single oscillator case. Let a' = r'e 10 , (3' = p'e 1 ^ , etc. Then, 

, 1/i; , l2 , „, 2 ,„ 2 , , l2 ,(r' 2 r" 2 + p ,2 p" 2 + 2r'r"p'p"cos(6)) m 

e -Kl-'T+lp'T+k'l 2 +|p'l 2 ) V UL IT \JJ- (5.40) 

where 6 1 = 9' — 9" — <p' + 0". The peak of this function is at 9 = 0, r' = r", p' = p", and 
r' 2 + p' 2 = m. Note that 9 = implies that the states evolve together as they do in 
the classical equation of motion ( |5.3| ). The correlation functions of the the position 
and momentum are give by 

WWQMcr) = + «") Walter), (5.4i) 



ma' 



(e/,/?"|Pil<W> = zVm, 1 ( W + / g/^/) - (^'^"l^')- ( 5 - 42 ) 

Likewise for the position and momentum of the second oscillator. Expanding about 
the peak, we have 

(a",/3"|4|0a'^> = -cos(0' - 0") + O(n) (5.43) 

(a",/?"|A|0a'/3') = r'sm(0' - 0") + (5.44) 

(a",P"\Q 2 \^r) = -cos(0'-0") + 0(fr) (5.45) 

a; 

(e/,/3"|A|<^> = p'sin(0' - 0") + 0(h) (5.46) 
The energy for each of the oscillators in the classical limit is 

{<t>a>p\El\<l>a>P>) = r' 2 + 0(H), {<j>a>P>\E2\<l>a>p) = / + 0(K), 

(<t>a'r\E t otal\<f>*'P>) = m ' = r ' 2 + P' 2 + 0{h) (5.47) 



So, we see that we meet the constraint equation (|5.5f) , and we get back the classical 



equations of motion (|5.3j) in the classical limit 
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6 Discussion 



It is interesting to note that the gauge degree of freedom for both of these sys- 
tem comes out in a phase factor in front of the reduced phase space coherent state. 
Although is might just be a artifact of the harmonic oscillator (s), it is suggestive 
that the gauge should appear in this way in general. If this is the case, the only 
dependence on the gauge orbits in the path integral will appear in the one form, 



(p>Q\i\P>q)dt = pdq 

(p, q\ e - inp ' q) ^e lfM \p, q)dt = pdq + df (6.1) 
at 

If there is not a boundary, this difference can just be integrated out as a total deriva- 
tive. If there is a boundary, we pick up a boundary term that is still dependent on 
the gauge orbits. For example, such a gauge symmetry break term is seen in the 
relationship between the Chern-Simons actions and the Wess-Zumino-Witten action. 
It should be noted that, the above correlations functions ( 4.18|) and (|5.41|) are 



not physical observables. This correlations are still dependent on the gauge degree 
of freedom. However, they do show that the equations of motions are still embedded 
in the formulation of the coherent state on the reduced phase space. In terms of the 
double harmonics oscillator, it is interesting to note that the width of the correlation 
function in the angular direction is dependent on the energy of each oscillators, 

E total <-.\ 

In terms of "quantum clocks," this means that at low energies, the correlations be- 
tween this clocks may become fuzzy. Certainly, a more precise statement in terms of 
observables needs to be considered. 
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